Formation of three-particle clusters in hetero-junctions and MOSFET structures 
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A novel interaction mechanism in MOSFET structures and GaAs/AlGaAs hetero-junctions be- 
tween the zone electrons of the two-dimensional (2D) gas and the charged traps on the insulator 
side is considered. By applying a canonical transformation, off-diagonal terms in the Hamiltonian 
due to the trapped level subsystem are excluded. This yields an effective three-particle attrac- 
tive interaction as well as a pairing interaction inside the 2D electronic band. A type of Bethe- 
Goldstone equation for three particles is studied to clarify the character of the binding and the 
energy of the three-particle bound states. The results are used to offer a possible explanation of the 
Metal-Insulator transition recently observed in MOSFET and hetero-junctions. 
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Significant advances in low temperature physics are 
particularly connected with the recent successes in fabri- 
cating two-dimensional (2D) electronic structures char- 
acterized by high mobility. The formation of the inver- 
sion layer of particles on interfaces allows the study of 
the unusual low-temperature behavior of a 2D electron 
gas , as well as the examination of theoretical predic- 
tions. Particularly, experimental measurements at very 
low temperatures, T ^ 2K, show quite nonstandard re- 
sults, like the fractional quantum Hall effect^ and the 
Metal-Insulator transition (MIT)^'^, the origin of which 
is not yet properly understood. It is worthwhile to ob- 
serve that the fractional quantum Hall effect and MIT 
occur only under some special regimes where the temper- 
ature and the impurity concentration in the samples are 
very small. On the other hand the observation of such 
unusual effects only in GaAs/AlGaAs hetero-junctions 
and Si-MOSFET raises the question whether structural 
peculiarities of these devices are responsible, or rather 
some fundamental law is observed. 

The band structures of GaAs/AlGaAs hetero- 
junctions and of Si-MOSFET are well known, and 
their general features are identical despite their differ- 
ent structures. The inversion of current carriers oc- 
curs at the interface of two semiconductors (or insu- 
lator/semiconductor) with different band gaps. These 
semiconductors are exclusively doped by p- and n-type 
impurities in order to get high mobility in the samples. 
In extremely clean samples, the doping by one type of 
either acceptors or donors forms single levels in the gap 
of each semiconductor. In the process of, e.g., electron 
inversion, the donor centers in Si02 of the MOSFET 
structure (or in, Al^Gai^xAs hetero-j unction) become 
positively charged by transferring electrons to the 2D 
electronic band. The charged donors are located within 
a region of at most ~ 200A'* from the oxide-silicon in- 
terface, and their energy level lies above the Fermi level. 
Such interfacial charged states act as trap centers for the 
band electrons. The density of trapped centers in, e.g., 



MOSFET structures is of the order of lO^cm"^^, which 
is considerably smaller than the typical carrier concen- 
tration of ~ 2 X lO^^cm^^. The Coulomb potential of 
charged traps seems to be screened imperfectly due to 
their low densities, and the scattering of the band par- 
ticles on these trap centers may be more essential than 
the intra-band particle-particle scattering. In this letter 
we study the effects of the scattering of band electrons 
on the charged traps. 

The Hamiltonian of the model can be written in the 
form H = Hq + Hint , where 
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In Eqs.(l) and (2), ^ and b'^^ (ak^o- and 60,0-) are the 
creation (annihilation) operators for the electrons in the 
band and in the trapped levels, respectively. 

The trap centers are modeled for simplicity as a dis- 
persionless single level. e(k) and wq in Eq.(l) are the 
energies of the band electrons and of the trap level, re- 
spectively. The energy of the trap centers is considered 
to be larger than the chemical potential fj, of the band 
electrons, ujq > fj.. The first term in the Hamiltonian 
Hint represents the scattering of two band electrons via 
the interaction potential ^(k) followed by the trapping 
of one of them by the donor level. The second term rep- 
resents the scattering of a band electron with a trapped 
electron (on a donor level) with finally turning both of 
them into the band. 

For temperatures kT < {ujq — fi), the trapping cen- 
ters contain a definite number of electrons at thermody- 
namic equilibrium, and the considered trapping mecha- 
nism is assumed to be essential. However, the mechanism 
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is destroyed by increasing the gate potential, for then the 
chemical potential /i reaches the trap level, causing the 
filling of all trapping centers. A temperature increase 
also leads to the destruction of the trapping mechanism. 

We apply a unitary transformation^ with the inten- 
tion to result in the cancellation of the off-diagonal term, 
given by Eq.(2). We will show here that this unitary 
transformation creates a "trap mediated" effective at- 
traction between three electrons in the band^. Notice 
that our approach to the problem is similar to the can- 
cellation of the phonon subsystem in superconductivity, 
where a canonical transformation yields an effective at- 
traction between electrons, see, e.g.,^. Let us expand a 
new Hamiltonian H = e~^^He^^ in power series of the 
operator S* = 5*1 + 5*2 -I- 5*3 -I- • • •. Regrouping the terms 
of the same order in Hint, the conditions which define 
Si, i = 1,2, ■■ ■ are obtained recursively. Si is determined 
by Hint — i[Si, Ho] = and leads to 
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The new Hamiltonian H can now be written in the form 
H = Ho — ^[Si,Hint] — i[S2,Ho]. S2 is obtained from 
the condition that the off-diagonal terms in the Hamil- 
tonian be cancelled. The equation for H then becomes 

H = Ho — ^[Si, Hint]diag, where the last term contains 
only diagonal ck;nicnts. Introducing Eq. (3) into the ex- 
pression for H finally yields the effective Hamiltonian 



H = Hq + He-e + 
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where Hq, Hp-e and i^A describe the effective one- 
particle Hamiltonian, the electron-electron interaction, 
and the three-particle clustering, respectively. The one- 
particle effective Hamiltonian Hq is given by 

Ho = Ho -'^b+^X,,7''^£iiqL)a+^aq,a 
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In thermodynamical equilibrium < b^^,bo,a >= Scr,a', 
and the one-particle energy is renormalized, Hq = 
E ^(q)a;i',^aq,^ with e(q) = e(q) + ei(q) + J(q). 

The electron electron interaction Hamiltonian He-e 
also contains terms with spin flipping due to the exchange 
scattering of 2D electrons on trapped ones. In thermo- 
dynamic equilibrium, -ffg-e has the usual form 
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where the effective two-particle interaction potential 
K-e'^' (ki , k2 , q) appears to be attractive. 
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Indeed, the denominator of each term in v}l{.^^ is nega- 
tive, owing to the fact that a donor level lies higher than 
the chemical potential of band electrons. 

The third term in the Hamiltonian (4) describes an 
effective three-particle scattering, 
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Using again the condition of < wq, it is possible to 
see that the strength of the three-particle interaction is 
negative, which results in formation of clusters of three 
electrons. 

This effective attraction among three electrons can be 
understood according to the following physical argument. 
The proposed mechanism of two-particle interaction with 
trapping of one of the particles, in contrast to an intra- 
band electron- electron scattering, destroys locally the 
electro-neutrality of the 2D electron gas. The necessary 
electro-neutrality in hetero-junctions or in MOSFET's is 
restored by the ensuing adaption of the height of the 
Schottky barrier, i.e. by a change in the value of the 
band bending energy. However, the trapping and releas- 
ing processes are so fast that the barrier's height cannot 
follow. As a result of the trapping of band electrons, a 
hole appears which acts as an attractive center for other 
electrons. 

The energy level of the trap centers in the above calcu- 
lation is chosen to be dispcrsionlcss for simplicity. How- 
ever, even in the single level case, the donor center ener- 
gies depend on the spatial coordinates of the impurities 
due to the band bending, and therefore become disper- 
sive. Including the dispersion of the trap level does not 
change ciualitatively our results. 

We now proceed to show that a three-particle attrac- 
tive interaction can lead to the formation of a bound 
state. To this end, we consider for simplicity only the 
three-particle interaction, and neglect the pairing in- 
teraction. The Schrodinger equation for three identi- 
cal particles in 2D with a generic interaction potential 
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of the form V{ri — r2; ra — ri; r2 — ra) is written in the 
foUowing form after introducing the Jacobi coordinates 
R = i(ri + r2 + rs), r = rg - ri and z = 



- r2 



xV'(R,r,z) = (e + 3eF)V(R,r,z), (9) 

where e is the three-particle excitation energy measured 
from the three-particle Fermi level. After excluding the 
center of mass coordinate R by expanding ^/'(R, r, z) in 
plane waves, V(R,r,z) = EQ,p,^ei^'^^+P''+^'^^{p,q) , 
a Bethe-Goldstone-type equation, similar to the equation 
for Cooper pairs^, is obtained. 
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The interaction potential V{p\ q'; p, q) is assumed to be 

attractive when the energies of the three particles (before 
the coordinate transformation) lie in a narrow vicinity 
fiu>o of the Fermi surface {cft^f + Tiuiq). Here Twoq is a 
cut off energy which is comparable to the order of trap 
level energy measured from the Fermi level. This condi- 
tion restricts the energies, jep < /2m < |(eF -|- Twoq) 
and €p < (f' jlm < e^? + hujQ. of the quasiparticlcs ob- 
tained after coordinates transformation. Therefore, the 
simplified attractive interaction for a system of linear 

size L is y(p',q';p,q) = for fe^ < C < 

|(eF + ?iwo), < ■^■,^<iF ^ huo, and otherwise. 

The following equation for the bound state energy of 
a three-particle cluster results, 
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where is the value of the density of electronic states 
on the Fermi surface. The integration of Eq.(ll) gives 
an equation for e, 

^ o, 6-2? -4e(3-e) 
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where e = e/fkuQ is the dimcnsionlcss excitation energy. 
The study of Eq. (12) for arbitrary negative values of 
e shows that the r.h.s. of this equation is a monotonic 
and positive function with a maximum value equal to 
In 8 at e = 0. This imposes a lower restriction on the 
attractive potential, Vq > in2Jv^yia>o " Consequently, for 
attractive potentials strong enough, Eq.(12) possesses a 
unique negative solution for — e € (— oo,0). This implies 
the existence of three-particle bound states. 

In a weak coupling regime, when — e <C 1, Eq. (12) is 
simplified to 
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the solution of which does not possess a gap-like form 
for the excitation energy. 

In the opposite case of large negative solutions,—? > 3, 
Eq. (12) leads immediately to the following result for 
the bound energy in the strong coupling regime, e = 
— jVoihuioN)"^ , which shows a clear perturbative non- 
collective behavior. 

This discussion shows that molecular clustering rather 
than a coherent state is realized in the system. The 
ground state of the system becomes unstable with re- 
spect to the three-particle attraction. This seems to lead 
to molecular type formation with negative energy. 

In conclusion, we want to emphasize that it is possi- 
ble to understand qualitatively the reason of MIT occur- 
ring at very low temperatures in GaAs/AlGaAs hetero- 
junctions and MOSFET^'^ in the framework of the for- 
mation of three-particle bound states we describe above. 
The elastic scattering of electrons on impurities at low 
temperatures, which is characterized by a relaxation time 
To, results in the localization of all electronic statcs^°, un- 
der the condition that h(jJo{n) < kT < h/ro, producing 
an insulating behavior for conductivity. Here uioin) is the 
trap level energy measiircd from the Fermi energy, which 
is a function of the 2D electron concentration n, or the 
gate potential. On the other hand, in a regime corre- 
sponding to kT < hLOo{n) < h/ro, which can be reached 
by varying the electron concentration or the tempera- 
ture, the formation of three-particles bound states results 
in the vanishing of the weak localization corrections to 
conductivity. This is due to the fact that the scattering 
of the three-particle clusters on the impurities does not 
lead to quantum interference. Instead, the cluster's wave 
function accumulates an additional phase by rotation of 
the cluster in the process of scattering, while the cen- 
ter of mass motion of the cluster is still extended. The 
expression for the conductivity can be written as 
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where the first and the second terms in Eq. (14) 
correspond to the Drude and the weak localization 
contributions,^^, which correspond to three-particles 
clusters and free band electrons with concentrations of 
riA and n/ respectively, is the inelastic scattering 
time Tin = aT^P where a is some constant,^ > 2 and 
p = 2 for probable electron-electron scattering mecha- 
nism. Notice that a logarithmic temperature dependence 
of cr in the metallic phase has been observed,^, in high- 
mobility n — Si — MOSFET which is in good agreement 
with our assumption. Observation of the negative low- 
field magnctorcsistance in the metallic phase also sup- 
ports an important role of the quantum interference ef- 
fects in MIT. We neglect in Eq.(14) an additional loga- 
rithmic quantum correction due to the electron-electron 
interactions,^^, which is responsible for positive magne- 
toresistance also observed in experiments, 2. 
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Charge conservation allows us to write the total con- 
centration of particles as n = + n/ + nj, where 
nt is the concentration of trapped electrons, which 
exponentially decreases with temperature as nt(T) = 
rij exp{— a;(n)/A:T} with n° being the concentration of 
trapped impurities. Assuming that the clustering oc- 
curs at T = Tcin), ua can be expressed near Tc{n) as 
^ = We rewrite Eq. (14) in the form 

^ = l + 2^-^ln^, (15) 

where (tq = e^Ton/m is Drude conductivity, and T* = 
y/a/TQ. In Eq.(15) wc neglected the trap level contribu- 
tion due to nt < n. Rescaling T by T* as T/T* = t 
for r < 1 and choosing the parameter of randomness 
A = 2iTepTa ^ ^ ' unknowing parameter Tc/T* = Tc 
in Eq.(15) can be extracted by fitting a{T)/ao to the ex- 
perimental data. The temperature dependence of cr is 
drawn in Fig.l for the best fit parameters Tc{n). 
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FIG. 1. The temperature dependence of the ex- 
perimental (dots) conductivity scaled to the critical 
oiic^^. The densities increase from below to above as 
n = 6.85, 7.17,7.25, 7.57, 7.85x10"' cm"^ The lines are the- 
oretical fits according to (15) where the lower three curves 
are assumed to be undercritical (without tia) with the fit 
Tc = 3.9,3.6, 1.7K and A = 0.330,0.078,0.004 correspond- 
ingly and the above two metallic curves arc critical (with ua) 
with the fit Tc = 1.4, 0.8K and A = 0.48, 0.86. 

While the overall fit to MIT is satisfactorily we see 
deviations on the insulating side. This seems to be due 
to strong localization effects in these experiments. Ef- 
fects observed in the weak localization regime,^, showing 
a more pronounced logarithmic behaviour. The value of 
the conductivity at the critical density Gc observed in 
the two experimets strongly differs from each other: in 
the strongly disordered case,^, Gc ~ e^//i, whereas in 
the weakly disordered case of"^ Gc ^ 120e^/h. Therefore, 
we think that the conductivity at the critical density Gc 
does not show a universal behavior. It depends on two 
factors: on the impurity concentration in the Si substrate 
or in the 2D electron gas and on the electron concentra- 
tion (or on the band filling) in the inversion layer. By 



increasing the band filling in the insulator side of the 
MIT the three particle clusters appear which weaken the 
localization tendency in the 2D electronic system since 
the scattering off the cluster on the impurities does not 
lead to the quantum interference effects. At the criti- 
cal density, the contributions coming from the clustering 
completely compensate the localization corrections and 
the conductivity is defined by the value of the residual 
Drude conductivity, which is temperature independent at 
low temperature. 

Notice that another mechanism for MIT, which is 
controlled by a temperature-dependent trapped-electron 
concentration ni{T), has been recently proposed by Alt- 
shuler and Maslov^"^. As the comment and reply shows, 
this qualitatively correct explanation cannot reproduce 
quantitative features of the experiment. A critical dis- 
cussion of different approaches can be found in Ref.^. 

Although we have not discussed a role of effective pair- 
ing governed by Eq. (6) in the Hamiltonian, there was an 
attempt to interprcte the experimental data on MIT as 
a result of possible superconducting ground state, It 
is well known that the effective pairing is suppressed by 
the order parameter phase fiuctuations in 2D systems 
reducing Tsc of the superconducting transition to zero. 
However, fluctuations of the order parameter modulus 
above Tsc may lead to the metallic phase. 

The geometry of the cluster may be either in the form 
of a triangle with 3/2 and 1/2 total spin, or of a string- 
like configuration with 1/2 total spin, when two electrons 
with antiparallel spins are placed at the same point and 
the third electron with arbitrary spin is far from them. 
In the case of 3/2 total spin, a magnetic field parallel to 
the triangle area does not destroy the cluster, whereas 
the configurations with 1/2 total spin are destroyed due 
to the Zeeman effect. In both cases the magnetic field 
effects are defined by the contributions coming from the 
quantum localization corrections. 

The model of three-particle clustering due to the dis- 
cussed exchange type of interaction with donor levels 
seems to be also a favorable candidate for the understand- 
ing of the Fractional Quantum Hall effect. The interac- 
tion of the band electrons with trap centers effectively 
leads to a formation of three-particle clusters, see Eq. 
(8), as well as to the superconducting fluctuations due to 
effective pairing interactions, Eq. (6). Both mechanisms 
decrease the ground state energy of the system. Strong 
magnetic fields in the quantum Hall regime polarize the 
spins of molecular clusters and a triangular geometry for 
the cluster is realized due to the Pauli principle. An anti- 
symmetric orbital wave function of the triangular cluster 
will contain ab initio the Jastrow prefactor. The angular 
momentum M = 3 of the cluster provides a natural ar- 
gument in Laughlin's theory to connect the filling factor 
u = 1/3 of the parent states with the angular momentum 
M = 3. 
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